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Summary 


Presents  a  linear  solution  for  determining  the  effect  of  combined 
axial  and  lateral  loads  under  which  a  cylindrical  sandwich  shell 
will  buckle.  The  facings  of  the  sandwich  cylinder  are  treated  as 
homogeneous  isotropic  cylindrical  shells  and  the  core  as  an  ortho¬ 
tropic  elastic  body.  The  characteristic  determinant  that  represents 
the  solution  to  the  problem  is  solved  numerically.  Curves  are  given 
that  show  how  the  buckling  load  changes  as  the  parameters  of  the 
problem  change. 
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Introduction 


Sandwich  construction  is  a  result  of  the  search  for  a  strong,  stiff,  and  yet  light  weight 
material.  It  is  usually  made  by  gluing  relatively  thin  sheets  of  a  strong  material  to  the 
faces  of  relatively  thick  but  light  weight,  and  often  weak,  material.  The  outer  sheets  are 
called  facings  and  the  inner  layer  is  called  the  core. 

Such  a  layered  system  presents  difficult  design  problems.  What  is  offered  here  is  a 
straightforward  method  for  dealing  with  some  of  these  problems. 

The  problem  to  which  the  method  is  applied  is  that  of  the  elastic  stability  of  a  sandwich 
cylinder  under  uniform  external  lateral  load  and  uniform  axial  load. 


Notation 


r,  6,  z  radial,  tangential,  and  longitudinal  coordinates,  respectively 

a  radius  to  middle  surface  of  outer  facing 

b  radius  to  middle  surface  of  inner  facing 

r  *  mean  radius 

m 

t  thickness  of  each  facing 

—This  ResOarch  Note  is  a  revision,  under  the  same  title,  of  Forest  Products  Laboratory 
Report  1852,  issued  in  1955.  It  was  originally  prepared  by  Everett  E.  Haft,  and 
issued  as  one  of  a  series  by  the  Forest  Products  Laboratory  i  r,  cooperation  with  the 
U.S.  Navy,  Bureau  of  AeronauTics. 
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■^Maintained  at  Madison,  Wis.,  in  cooperation  with  the  Universily  of  Wisconsin. 
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length  of  cylinder 
modulus  of  elasticity  of  facings 
Poisson’s  ratio  of  facings 
modulus  of  rigidity  of  facings 

modulus  of  elasticity  of  core  in  direction  normal  to  facings 
modulus  of  rigidity  of  core  in  r@  plane 
modulus  of  rigidity  of  core  in  rz  plane 
intensity  of  uniform  external  lateral  loading 

I 


I  +  >  -  In  * 

a  E  a  a 
c 


Tr0’  Trz 


u.  V,  w 


m 


r0 


rz 


normal  stress  in  core  in  radial  direction 
transverse  shear  stresses  in  core 

radial,  tangential,  and  longitudinal  displacements,  respectively 
number  of  waves  in  circumference  of  buckled  cylinder 
number  of  half  waves  in  length  of  buckled  cylinder 
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N  ,  N  ,  N  normal  forces  and  shear  force  per  unit  length  of  facing 
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transverse  shear  forces  per  unit  length  of  facing 
bending  moments  per  unit  length  of  facing 
twisting  moments  per  unit  length  of  facing 
surface  forces  per  unit  area  of  facing 
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(a  -  b) 


the  value  of  the  lateral  load  parameter  4^  when  the  axial  load  parameter 
<f>2  is  equivalent  to  zero 

the  value  of  the  axial  load  parameter  <j>  when  the  lateral  load  para¬ 
meter  <f>^  is  equivalent  to  zero 

natural  logarithm 


2r 


b 

a  or  1 1  +  A 


total  axial  load 


A,  B,  C,  D,  K,  L,  A' ,  B\  A",  B"  arbitrary  constants 

note  —  any  of  the  above  terms  that  appear  with  a  prime  (as  N  1 )  refer  to  the  inner 
facing. 
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OUTER  FACING 


Z  M  105  662 


Figure  I. — Sandwich  cylinder. 


a  im  «n 


Figure  2. — Differential  elements  of  core  end  facings  before  buckling 
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Mathematical  Analysis 


As  previously  stated,  the  core  is  relatively  weak.  Because  of  the  high  strength  of  the 
facings  the  core  need  carry  little  tension  or  compression  except  in  a  direction  perpen¬ 
dicular  to  the  facings.  The  facings  are  able  to  resist  shearing  deformation  in  their  plane 
and  it  is  necessary  only  that  the  core  be  able  to  resist  shear  in  the  radial  direction  in 
planes  perpendicular  to  the  facings.  In  this  analysis,  the  core  is  considered  to  be  an 
orthotropic  elastic  body.  It  is  unable  to  resist  deformations  other  than  those  just  men¬ 
tioned. 

Interdependence  of  the  core  and  the  facings  is  gained  by  equating  their  displacements 
at  the  interfaces.  To  simplify  the  analysis,  the  core  is  assumed  to  extend  to  the  middle 
surface  of  each  facing. 

Figure  1  shows  the  cylinder  and  the  coordinates  that  are  used. 


Prebuckling  Stresses 

Before  buckling  occurs,  the  cylinder  is  in  a  state  of  uniform  compression.  The  axial 

load  is  carried  by  the  facings  since  the  core  material  is  assumed  to  be  incapable  of 

carrying  load  in  this  direction.  With  facings  of  like  material,  the  stress  is  the  same  in 

both  facings.  Since  the  facings  have  the  same  thickness  then  the  loading  per  unit  length 

of  facing,  N  or  N  ' ,  will  be  the  same.  This  means  that  for  a  total  load  P 
z  z  — 

27iaN  +  2itbN'  =  P 
z  z 

The  calculation  of  stresses  due  to  the  lateral  pressure  is  a  problem  in  rotational 
symmetry.  Differential  elements  of  the  core  and  of  the  facings  are  shown  in  figure  2. 

Summing  forces  in  the  radial  direction  gives  for  the  core 


for  the  outer  facing 

aq  -  a  (CM  -  N0  =  0, 

r  =  a 

and  for  the  inner  facing 

b(o  )  -  N1  =  0. 

r  r  =  -  6 

Since  o  =  E 
r  c 
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Nfi  =  f.t  (  -)  ,  and 

v  3 

r  =  a 


mi  -  p*  (  — i  .  these  equations  can  be  solved  for  o  ,  NQ,  and  N ' 

N0  '  a  r  =  b  _£  _JL  -J 

3 

The  results  are- 


N’ 

b 


q?k 

qa( I  -  k) ,  and 
qak  where 


h  Et  ,n  T 

+  £  -  — = - - 

a  E  a 


As  P  and  q  Increase,  N  ,  N',  N  ,  N  ’  and  o  also  increase.  Eventually  a  condition  may 
*“•  _ z_  z  v  U  r 

be  reached  where  a  slight  increase  in  load  causes  the  cylinder  to  lose  its  state  of  uniform 
compression  and  buckle  as  a  result  of  elastic  instability.  This  buckling  is  assumed  to 
cause  only  a  small  change  in  the  stress  distribution.  These  small  changes  will  now  be 
considered. 


Buckling  Stresses 

The  core. — A  free  body  diagram  of  an  element  of  the  core  is  shown  in  figure  3. 

Neglecting  terms  which  are  products  of  more  than  three  differentials,  a  summation 
of  forces  in  the  r,  J  and  z  direction  gives 


3or  „  3t 

a  +  r  ~  +  +  ~SL  =  0 

r  3r  30  3z 


3t 


re 


„  +  2T  =  0 
3r  lrQ 


3t 


T  +  r  - 
rz  3r 


rz 


=  0 


Equation  (2)  may  be  integrated  to  give 


Ve'7  V9’  Vz’ 


(1) 

(2) 

(3) 

(4) 


■^Ravllle,  M.  E.  Analysis  of  long  cylinders  of  sandwich  construction  under  uniform 
external  lateral  pressure,  Forest  Products  Laboratory  Report  1844,  1954. 
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Figure  3. — Differential  element  of  deformed  core.  Z  M  103  364 


Equation  (3)  may  be  integrated  to  give 
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,  and  t  as  defined  in  terms  of  u,  v,  and  w  are 
re  rz  ~ 


a 

r 


=  E 


3u 
C  3r 


(5) 


(6) 


T 


re 


T 


rz 


rz 


3v 

3r 


3w 

3r 


v 

r 


(7) 

(8) 


It  is  convenient  to  assume  the  displacements  u,  v,  and  w  in  the  form 


u  =  f„(r)  cos  ne  cos-  z 

1  a 

(9) 

v  =  f  (r)  sin  n0  cos  -  z 

<L  3. 

(10) 

w  =  f  (r)  cos  ne  sin-  z 

O  Jl 

(11) 

This  form  will  permit  a  unique  determination  of  f  (r),  f  (r),  and  f  (r);  assumes  n  cir- 

X  £»  O 

cumferential  waves  and  m  longitudinal  half  waves  upon  buckling;  result?  in  zero  dis¬ 
placements  in  the  radial  and  circumferential  directions  at  the  ends;  and  imposes  no 
moment  upon  the  facings  at  the  ends. 
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Fiaure  4. — Differential  elements  of  deformed  facings. 
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The  facings. — Sinoe  the  problem  of  stability  of  homogeneous  cylindrical  shells  has 
been  solved  by  other  authors,  it  is  only  necessary  in  the  stability  analysis  of  cylindrical 
sandwich  shells  to  consider  it  as  a  composite  of  two  homogeneous  shells  bonded  together 
by  an  elastic  core  subject  to  the  compatibility  requirements  of  equal  displacements  at 
their  interfaces,  which,  to  simplify  the  analysis,  is  assumed  to  extend  to  the  middle 
surface  of  each  facing. 


Free  body  diagrams  of  a  facing  element  showing  the  sense  of  the  forces  and  moments 
are  found  in  figures  4  and  5. 


Moz 


Z  M  106  V)2 

Figure  5. — Differential  element  of  facing  showing  moments  and  twists. 


The  equations  of  equilibrium  for  the  facings  can  be  obtained  from  a  mathematical 
theory  of  thin  shells.  It  is  felt  that  for  a  small  deflection  analysis  the  theories  presented 
by  various  authors  differ  principally  in  second-order  effects,  as  is  evidenced  by  com¬ 
paring  the  work  of  Flugge^  and  Timoshenko^  on  the  buckling  of  cylindrical  shells. 


a 

-flUgge,  W.  Stresses  in  Shells,  Springer-Verlag,  Berlin,  1962. 

5 

—Timoshenko,  S.  Theory  of  Elastic  Stability,  McGraw-Hill,  1936. 
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As  a  result,  the  equations  of  equilibrium  as  presented  by  Timoshenko  and  subject  to 
the  necessary  transformations  will  be  used  here.  Such  theory,  when  applied  to  cylindrical 
shells  of  radius  r,  yields  the  following  differential  equations: 

3N  3?u  82v  »  / 3?u  _  3v\ 


o  ^  azse  dzj 
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302 

(25) 


+  If  -  %  <26) 


•  +  rO^  (27) 

i 

1  a2u  i 

r  302  rde) 

(28) 
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Assuming  that  resultant  forces  other  than and  are  small  and  neglecting  products 

of  these  forces  and  the  displacements  u,  v,  and  w  which  are  also  small,  and  assuming 
further  a  membrane  analysis  where  the  bending  moments  and  transverse  shear  forceB 

i  n  i 


in  the  individual  facings  are  neglected, 
reduced  to: 
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the  equilibrium  equations  can  be 

(29) 

(30) 

(31) 


Evaluating  equations  (29),  (30),  and  (31)  for  each  facing  or  when  pjL  and  Zl^nd 
satisfying  the  compatibility  requirements  at  these  interfaces  provide  the  six  equations 
necessary  for  the  determination  of  the  six  arbitrary  constants  found  in  expressions  20. 

21,  and  22. 


qa(l  -  k) 
res- 


Mathematical  Analysis 

N  and  N„  of  equations  (29).  (30),  and(31)  are  replaced  by  \  and  jq 

z  u  L _  m  J  —A 

+  N  1  ,  respectively,  for  the  outer  facing  and  by  —  +  N^j  and  j^qak  +  N0J 

nectively,  for  the  inner  facing  (primes  denote  inner  facing).  This  necessary  because 
ttie  forces  in  the  buckled  shell  are  the  prebuckling  forces  plus  the  forces  due  to  buckling. 
The  N  ,  N  ',  Nq,  and  N  ’are  small  forces  due  to  buckling  which  are  later  to  be  expressed 

in  terms~of ""displacements,  (which  are  assumed  to  have  differed  negligibly  from  the 
prebuckled  state)!  and  the  product  of  these  small  forces  and  the  displacements  are  also 
neglected.  Expressions  (29),  (30),  and  (31),  as  evaluatedfor  the  outer  facing  then  become. 
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(32) 


(33) 


(34) 
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As  is  customary  in  such  stability  problems,  the  stretching  of  the  middle  surface  is 
taken  into  account  by  assuming  that  the  prebuckling  stress  resultants  are  very  large  in 
comparison  with  the  other  stress  resultants;  thus  the  terms 


—  (I  +  ea) 
4irr  9 
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qa( I  -  k) ( ‘  +  e  ) 


should  be  substituted  for  the  quantities  - - andqa(l  -  k),  respectively,  and  the  surface 

4ttt 

forces  X,  Y,  and  Z  should  be  multiplied  by  (1  +  £^(1  +  e  )  in  equations  (32)  to  (34).  In 
these  expressions 

e  =  1  *1  +  “ 

9  a  39  a 

e  =  (35) 

z  9z 

Finally,  the  expressions  for  the  forces,  moments,  and  twists  along  with  the  relation¬ 
ships  for  the  surface  faces: 

?=-(Vz> 


Y  =  "‘Ve’ 


X  =  q  -  (q  -  k  +  or> 


are  expressed  in  terms  of  u,  v,  and  w  and  substituted  in  (32),  (33),  and  (34).  They  lead  to 
the  following  equations  for  the  outer  facing: 

_  A  _  2  32w  ,  (I  +  y)  _  a2v  .  .  9u  .  (I  -  U)  32w 
EFz  -  0  -  a  +  — —  -  + —JL. — 


--♦i(l  -k)  fzie  +  fz 


.2  „  \  2,  ,  2. 
9  v  9u  )  a  (  I  -  u  / 


G  |!U  !* 
rz  9z  9r 


—  -  «  -  Cl  +  u)  .  92w  (I  -  u)  2  92v  9u  2  92v  92v 

EF0  -  ° - 2  3  9?96  +  - 2  3  „  2  +  96  3  *2  .2  +  ..2 

oZ.  o <■  oo 

a2( I  -  y2)  r  / 1  9u  9v  v  \ 

Et  r0  (a  96  9r  a  I 

(2  v  2 

,  9  u  l  2,  9  *i  .  3u 

U  +  ^2  r  3  *2  ~2  -  38  9? 

where  and  $  are  as  defined  under  notations. 
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To  achieve  proper  interaction  between  the  core  and  the  facings,  the  displacements  of 
the  middle  surfaces  of  the  facing  are  set  equal  to  the  displacements  of  the  core  at  r  -  a 
and  r=b. 


Thus,  displacements  u,  v,  and  w  in  equations  (37),  (38),  and  (39)  are  replaced  by  equa¬ 
tions  (20),  (21),  and  (22)  evaluated  at  r  =  a.  In  this  manner,  three  equations  in  six  arbitrary 
constants  (A,  B,  C,  D,  L,  and  F)  are  written  for  the  outer  facing.  In  a  like  fashion,  three 
equations  are  written  for  the  inner  facing.  These  equations  can  then  be  written  in  matrix 
form  as: 


ail 

ai2 

a!3 

ai4 

ai5 

ai6 

a 

a 

a 

a 

a 

a 

21 

22 

23 

24 

25 

26 

a 

a 

a 

a  . 

a 

a„„ 

31 

32 

33 

34 

35 

36 

a 

a. 

a 

a 

a 

a,. 

41 

42 

43 

44 

45 

46 

a 

a 

a 

a 

a 

a 

51 

52 

53 

54 

55 

56 

a  _ 

a 

a 

a 

a 

a„„ 

61 
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66 

wher 

e  the 

a  *s 

_ii 

are  defined  in 

terms 

=  <  0 


(40) 


Equations  (40)  are  satisfied  if  the  constants  A ,  B ,  C,  E^  F,  and  L  are  all  equal  to  zero. 
This  represents  the  uniformly  compressed  circular  form  of  equilibrium  of  the  cylinder. 
A  buckled  form  of  equilibrium  is  possible  only  if  equations  (40)  yield  nonzero  solutions 
for  the  constants.  This  requires  that  the  determinant  of  the  coefficients  of  these  constants 
be  equal  to  zero. 


It  is  possible  to  find  simultaneous  minimum  values  (absolute  sense)  of  <J>^  and  4>2  for 

which  these  six  equations  will  be  satisfied  for  any  values  of  the  arbitrary  constants. 
Mathematically  this  means  that  for  such  a  combination  of  loads  and  material  parameters 
the  deflections  are  indeterminate.  The  shell  becomes  elastically  unstable  and  the  loads 
that  bring  about  this  condition  are  the  desired  critical  loads. 


Numerical  Computations 

A  literal  solution  of  the  sixth  order  determinant  for  the  eigenvalues  (<}>  ,<{>  )  is  not 

1 

feasible.  A  numerical  solution,  from  which  curves  may  be  drawn,  is  possible  if  a  digital 
computer  is  used. 

For  sandwich  construction,  particularly  for  that  utilizing  honeycomb  cores,  the  equa¬ 
tions  can  be  greatly  simplified  by  assuming  E  ~  °c.  This  substitution  is  common  with 

c 

this  type  of  construction  due  to  the  small  relative  displacement  occurring  between  the 
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sandwich  facings.  Indications  are  that  analyses  based  on  this  assumption  are  sufficiently 
accurate  in  most  ranges,  but  for  relatively  short  cylinders  having  weak  cores  such  an 
assumption  may  result  in  serious  error. 

Poisson’s  ratio  for  the  facings  can  be  taken  as  one-third  and  the  value  of  k  closely 
approximated  by  the  fraction  one-half  for  most  cylinders  where  the  radius  is  large 
compared  to  the  sandv/ich  thickness  h. 


Making  these  substitutions,  the  a^’s  in  equations  (40)  are  found  to  be 
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The  problem  is  now  reduced  to  developing  a  search  technique  whereby  minimum  values 

of  4> .  and  $  can  be  determined  for  given  values  of  V  „  ,  V  ,  A,  and  Y,  such  that  the 
1  2  ro  rz  _  ~ 

determinant  goes  to  zero. 


This  problem  was  programed  for  the  digital  computer  and  results  are  presented  in 
figures  6,  7,  and  8.  The  solutions  show  a  linear  relationship  between  <J>  and  4>  as  was 

5,6  —  — 

found  previously  for  isotropic  shells —  wherein  the  slopes  of  the  lines  are  dependent 
upon  values  of  number  of  buckles  producing  minimum  loads.  The  lowest  of  these  lines 
are  shown  in  figures  6,  7,  and  8  for  chosen  values  of  Y.  These  figures  were  then  used  as 
a  basis  for  plotting  figures  9,  10,  11,  13,  14,  16,  and  17,  which  clarify  the  effect  of  the 
interaction  of  the  loads  by  presenting  them  in  terms  of  ratios. 


The  effect  of  special  cases  such  as  the  hydrostatic  case  can  be  seen  most  clearly  on 
figures  6,  7,  and  8  by  the  use  of  a  load  ratio  line  where  the  slope  of  the  line  and  the 
desired  load  ratio  are  equivalent.  The  hydrostatic  case,  for  example,  has  a  load  ratio  of 
one-fourth  or 


Discussion  of  Results 

2  h2 

For  illustrative  purposes,  the  three  values  of  A  -  — ~of  0.0001,  0.00001,  and  0.000001 

4r 

were  chosen  to  represent  the  range  of  delta.  For  a  given  delta  then,  figures  are  presented 

which  show  the  effect  on  the  critical  loads  of  (1)  varying  the  value  of  Y  while  holding  the 

values  of  V  and  V  a  constant,  and  (2)  varying  the  value  of  V  and  V  while  holding 
rz  r9  rz  r6 

irr 

the  value  of  Y  =  —  a  constant. 

L 

6,7 

The  results  of  this  analysis  agree  favorably  with  those  previously  presented^-  for 

lateral  pressure  ard  axial  load  acting  separately.  These  values  then  represent  the  values 

used  for  and  <J>  in  the  figures. 

1  ** 


For  the  case  of  lateral  load  only,  it  was  found  that  the  expression  for  $>^  is  given  by 


+A2(n2-  I)  3+4 
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The  value  of  4>  «  can  be  determined  with  the  use  of  figure  16,  which  was  taken  from 

7  — 

previous  work,- 


— Kuenzi,  E.  W.,  Bohannar,  8.,  and  Stevens.  G.H.  Buckling  Coefficients  for  Sandwich 
Cylinders  of  Finite  Length  Under  Uniform  External  Lateral  Pressure.  U.S.  Forest 
Service  Research  Note  FPL-OICK.  1965.  Forest  Prcd.  Lat>. 

— Zahn,  John  W.,and  Kuenzi,  Edward  W .  Classical  Buckling  of  Cylinders  of  Sandwich 

Construction  In  Axial  Compression — Orthotropic  Cores.  U.S.  Forest  Service  Research 
Note  FPL-018.  1963.  Forest  Prod.  Lab. 
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Figure  7. --Effect  of  combined 
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Figure  9. — Effect  on  critical  loads  of  varying  parameter  Y. 


Figure  10, — Effect  on  crilical  loads  of  varying  parameter  Y. 


Figure  II. — Effect  on  critical  loads  of  varying  parameter  X 
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Figure  14. — Effect  on  critical  loads  of  varying  parameter  Y 
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Figure  16. — Effect  on  critical  loads  of  varying  parameter  Y» 
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